Detection of continuous variable entanglement without coherent local oscillators 
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We propose three criteria for identifying continuous variable entanglement between two many-particle sys- 
tems with no restrictions on the quantum state of the local oscillators used in the measurements. Mistakenly 
asserting a coherent state for the local oscillator can lead to incorrectly identifying the presence of entangle- 
ment. We demonstrate this in simulations with 100 particles, and also find that large number fluctuations do not 
prevent the observation of entanglement. Our results are important for quantum information experiments with 
realistic Bose-Einstein condensates or in optics with arbitrary photon states. 

PACS numbers: 03.67.Mn,03.75.Ge,03.65.Ud,42.50.Dv 



The study of the quantum properties of matter-waves is a 
rapidly developing field known as quantum atom optics QUI]- 
Already several experiments have observed non-classical ef- 
fects in ultra-cold gases, including the Hanbury Brown- Twiss 
effect for bosons yl, anti-bunching for fermions 0], sub- 
Poissonian number fluctuations [5], and density correlations 
from molecular dissociation JfJ and in the Mott-insulator 
regime in an optical lattice 0]. Although impressive achieve- 
ments, the experimental techniques utilized in these obser- 
vations are insufficient to detect quantum squeezing or en- 
tanglement. The demonstration of entanglement — which 
Schrodinger described as being the central mystery of quan- 
tum mechanics JH] — will be an important step towards quan- 
tum information applications of ultra-cold atomic systems. 

In quantum optical systems continuous variable (CV) en- 
tanglement can be demonstrated experimentally by measuring 
certain correlation functions of electromagnetic field quadra- 
tures, and finding that these violate an inequality for sepa- 
rability, [l(J, or an inequality | IT] for demonstrating the 
Einstein-Podolsky-Rosen (EPR) paradox 111 211 . For a number 
of systems, these quadratures can only be determined using 
homodyne or heterodyne measurement techniques, and re- 
quire a coherent state local oscillator — a highly occupied 
mode of the electromagnetic field that is a good approximation 
to the output of many lasers 111311 . Such measurements have 
led to the observation of optical squeezing lfl4l [l5ll and the 
EPR paradox with photons lfl6ll . and have been used to per- 
form quantum state tomography IU7I1 and continuous-variable 
teleportation 111 811 . 

In principle squeezed or entangled atomic fields can be gen- 
erated from atomic four-wave-mixing lfl9l l2fj[ 21], molec- 
ular disassociation 12211 . or by mapping photon statistics 
onto atoms 112311 . Entanglement generated in these situations 
can potentially be used as a resource for quantum informa- 
tion 12411 . However, to unequivocally demonstrate entangle- 
ment between matter-waves it will be necessary to perform 
measurements sensitive to the relative phase of atomic wave 
packets. In principle it is possible to measure matter-wave 
quadratures in direct analogy to the optical case using atomic 
measurements with a suitable local oscillator (phase refer- 



ence) 025112611 . of which the matter-wave equivalent is a Bose- 
Einstein condensate. However these are typically not large, 
with a maximum of 10 s particles 02711 . and the phase stability 
is compromised compared to a laser, as atomic interactions 
result in a mode shape and energy that depend on the particle 
number. Unfortunately, the usual correspondence between ho- 
modyne measurements and the field quadratures that is ubiq- 
uitous in quantum optics may be lost if the local oscillator is 
neither large nor coherent. 

In this paper we show how to detect CV entanglement with- 
out a coherent state as a phase reference. We introduce three 
new entanglement criteria based on homodyne measurement 
that require no assumptions about the quantum state of the 
local oscillator. These allow for the demonstration of insepa- 
rability and the EPR paradox with measurements that are fea- 
sible with current or foreseeable technologies with ultra-cold 
atoms. We give a numerical demonstration of the application 
of these criteria in an experiment with 100 particles. 

We begin by reviewing homodyne measurement methods. 
The two quadrature operators for the mode with annihilation 
operator d are canonically defined to be 



X,, 



and Y„ 



(1) 



These operators do not conserve boson number and are not 
detectable directly by intensity measurements. In optics, ho- 
modyne and heterodyne measurements allow access to these 
observables by using a beam splitter to interfere the signal 
with a local oscillator with well-defined phase and amplitude 
followed by intensity measurement(s) as illustrated in Fig. Q] 
Using the Heisenberg picture, if a = a in and b = b\ n describe 
the signal and local oscillator before the beam splitting, then 
output modes a out and 6 out are given by 



= ta ir 



and 6„„ t = 



t*b ir 



(2) 



where |i| 2 + |r| 2 = 1. Here we focus on balanced homodyning 
using a 50-50 beam splitter with t = r = 1/ V2. The final 
step of a balanced homodyne measurement is to measure the 
difference in the number of particles exiting the beam splitter 
ports d out and b out . The measured quadrature is rescaled by 



2 




System 2 




FIG. 1 : Homodyne measurements of the two systems. Signal modes 
1 and 2 are mixed with local oscillators; the resulting number differ- 
ences are a measurement of the quadratures X 1 and X 2 . 



It has been shown 12911 that for any local observables on the 
two systems, Aj and Bj, with commutator [Aj, Bjj = iCj, 
all separable states obey 

Var [i x +i 2 ] +Var[4+4] > {(C^ + \(C 2 }\, (6) 

where Var[A] is the variance of the observable corresponding 
to the operator A. For our system, it follows that the two 
systems are entangled when the inequality 



Varfi^ ± X 2 ] + VarfYj T %} > 



1 



(a t 1 a 1 )/(6 t 1 6 i ; 



+ 2 



1 



(a\a 2 ) / (b\b 2 



(7) 



the size of the local oscillator, according to 



X„ 



out^out ^out^out 



\ in m/ 



(6tS)i/2 ' 



(3) 



If the local oscillator is a coherent state \/3), where /3 is real 
and positive, and the coherent state is large (wb) = |/3| 2 ^> 
(a) a), then the difference between the two quadrature opera- 
tors X c and X m will be negligible. For a signal that is com- 
parable in size to a coherent local oscillator |/3| 2 ~ (a^a), 
the moments of the canonical quadrature observables can be 
determined from the measured observables using the result 



(XJ = (X C ), (Xl) = (X 2 c ) + (ala)/(tfb). 



(4) 



However, for local oscillators that are not coherent states the 
difference between X m and X c cannot be neglected. 

For matter-waves the equivalent of the beam splitter oper- 
ation is performed by Bragg, RF, or Raman pulses fl]] to in- 
terfere atoms with different momenta and/or internal states. 
However, the best BEC phase references will typically have 
less than 10 7 atoms. With the additional effects of atomic in- 
teractions and finite temperatures, it can be expected that typ- 
ical BECs in the lab will have greater phase and/or amplitude 
fluctuations than the coherent output of a laser. Therefore we 
must take full account of the difference between the canonical 
and measured quadrature variables when detecting entangle- 
ment in condensate experiments. 

Separability criterion: The observation of CV entangle- 
ment requires the simultaneous measurement of two quadra- 
tures, and so balanced homodyning requires two phase ref- 
erences. We follow the general method of Duan et al. |@], 
Simon ifioll and Hofmann and Takeuchi |28] to derive a sep- 
arability criterion for the two systems of interest, depicted in 
Fig-El System 1 consists of the signal a 1 and phase reference 
b\, and System 2 has signal a 2 and phase reference b 2 . Mea- 
surements are made of the (number-conserving) quadratures 



X; 



1/2 



Y; 



(5) 



with commutator 



-2iS j>k (blb 



j j 



is violated. For semi-classical local-oscillators the right hand 
side (RHS) of Eq. (Q approaches four, as originally derived 
by Duan et al. |0] and Simon [flOn . Note that this result is not 
the same as applying the coherent state corrections [Eq. ©] to 
the derivation of Duan et al. Jg}] and Simon ifioll — for further 
discussion see the Appendix. 

EPR criteria: The value of a quadrature of System 2 may 
be estimated from measurements of System 1 due to the exis- 
tence of correlations between them. For separable states, the 
accuracy of this estimation will be restricted by the Heisen- 
berg uncertainty limit. However, entangled systems may al- 
low one to make predictions seemingly better than this fun- 
damental limit, a fact related to the Einstein-Podolsky-Rosen 
paradox |H by Reid ifTTll . 

We denote the inferred value of X 2 by the measurement 
result of X 1 as X™ 1 ~ /(XJ, and similarly Y 2 lul = g(Yi), 
with the error of these estimates given by Aj 2 nf [X 2 ] = ((X 2 — 
Xf f) and Af n{ [Y 2 ] - ((Y 2 - f ) 2 ), respectively. The 
EPR paradox is demonstrated when the inferred predictions 
have better accuracy than allowed by the Heisenberg uncer- 
tainty principle 111 ill . This occurs when the inequality 



A 



A X 2. 



> 



1 - (a\a 2 ) / (b\b 2 ) 



(8) 



is violated. Although the optimal inferred quadrature val- 
ues could be in principle any function of the measured val- 
ues, for simplicity we restrict ourselves to linear functions, 
i.e. X' 2 nl = a + bX 1 and Y 2 ni = c + dY v The sim- 
plest guess can be motivated by the inseparability criterion 
above; if the variances in Eq. (|7]i were vanishingly small 
then we would infer that X 
Yt f = ±Y X T (Yi) 



inf 



TX X ± (Xj) + (X 2 ) and 
(Y 2 ). This results in inference errors 



Af nf [X 2 ] =Var[X 1 ±X 2 ] ! 
A 2 nf [F 2 ] =V ta [Y 1 TY 3 ]. 

Together with Eq. ([8]), this implies 



Var[X 1 ± X 2 ] + Var[yj T Y 2 ] > 2 



1 - 



(b\b 2 



(9) 



(10) 



This is the same result as for the separability criterion above, 
Eq. ( [Tol l, except the RHS is twice as small (typically). Thus 
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it is more difficult to demonstrate the EPR paradox than in- 
separability, but it proves the stronger result that either local 
causality is violated or quantum mechanics provides an in- 
complete description of System 2 (for a discussion of this hi- 
erarchy see HEl). The RHS of Eq. CO) becomes two for 
semi-classical local operators which agrees with the work of 
Reid et al. $3$\ . This inequality is a special case of a more 
general EPR criterion for arbitrary observables lO . 

It is possible to derive a second EPR criterion that is vio- 
lated for a greater range of states than Eq. ( TTOb . The optimal 
linear inference of X 2 based on the measurement of X l is 111 111 



Fock state |100) 



Thermal (N) = 100 



Var [A] 



A - (A: 



(id 



where the covariance Var [A, B\ = (AB) — (A)(B). In sep- 
arate measurements we can estimate Y 2 a f ter measuring Y x in 
the same way. The inference errors from this method are 

Af nf [X 2 ] = Var [X 2 ] - Var [X, , X 2 ] 2 /Var [Xj , 
&L[Y 2 ] = Var[y 2 ] -Var^^f/Var^], (12) 

which are less than or equal to the comparable results in 
Eqs. (|9]). Substituting Eqs. ( fl2l i directly into Eq. ^ results in 
an EPR criterion that is superior to Eq. (TTOb . This result agrees 
with that of Reid ll ill in the limit that the phase-reference is 
large and the RHS of the inequality is unity. 

Numerical demonstration: We illustrate the use of these 
criteria in degenerate four- wave mixing with atoms |T^, 20ll . 
where condensate atoms in mode collide and populate 
modes 1 and 2 respectively. The Hamiltonian is 



H = ix(aoa\a J 2 - a^a^), 



(13) 



where is the annihilation operator for mode j and \ repre- 
sents the strength of the coupling. This model can be realized 
with a moving condensate in an optical lattice, and the outgo- 
ing modes 1 and 2 are predicted to be entangled 12 111 . 

Firstly, we consider an initial state with TV bosons in mode 
and vacuum in 1 and 2 (i.e. |*(0)) = |7V) |0) |0>). The 
Hamiltonian then ensures the state can be written in the form 

\*{t)) = J2 c ™{t)\N -2m)\m)\m), (14) 
m 

at all times. The Hamiltonian in this basis can be numerically 
diagonalized (for N < 10 4 ), providing a solution for the exact 
dynamics of the system. Fig.[2|a) shows the number of atoms 
in each mode as a function of time for TV = 100. 

To measure the quadratures of the signal modes (1 and 2) 
we require a phase reference for each. We can use the remain- 
ing condensate in mode as our phase reference, but must 
first divide it in two with a beam-splitter according to 



2, 



(15) 



where the auxiliary mode 3 is initially in the vacuum state. 
Remember that these local oscillators are not coherent states. 
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FIG. 2: (Color online) The population dynamics and entanglement 
criteria for the case of (a-c) exactly 100 particles and (d-f) a thermal 
distribution with a mean of 100 particles. In (a,d) we see that parti- 
cles are transferred from the condensate to modes 1 and 2. In (b,c,e,f) 
the entanglement inequalities are violated when the solid line enters 
the shaded areas enclosed by the dashed lines. as described in the text. 
The horizontal dotted lines indicate the inequality violation boundary 
for the canonical quadratures (^. [TolfTlll . 



We then perform the quadrature measurements as expressed 
by Eq. (O, and compare the results with the separability cri- 
terion, Eq. (0, and the EPR criteria, Eq. ( fTOt and Eqs. ( 181121 ). 
The results are plotted in Figures |2jb,c) where we can see that 
all the inequalities are violated at some stage in this experi- 
ment. After some time the number of particles in the signal 
beams grow larger than the phase references, and this mea- 
surement scheme is no longer optimal to detect the entangle- 
ment. 

The weaker EPR criterion, Eq. ( fTOt in Fig. |2jb), shows a 
smaller region of \t where violation occurs than the stronger 
version, Eq. ( fT2l in Fig. |2|c), as expected. The results in 
Figs. |2|b,c) identify times where it is necessary to fully treat 
the quantum nature of the quadrature measurement. Perform- 
ing these measurements and then using the criteria derived 
by Duan et al. 0], Simon ifioll and Reid 111] would lead to 
falsely identifying entanglement when the size of the phase 
references and signal beams become comparable. 
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A particularly interesting result is that the inseparability cri- 
terion [Eq. (0] identifies entanglement in a small region [the 
top-right corner of Fig. |2](b)] when the local oscillator is ac- 
tually smaller than the signal modes and has sufficient fluctu- 
ations that the original Duan criterion ana would not detect 
entanglement. This effect becomes more prominent for larger 
numbers of particles (see the Appendix for further details). In 
this regime it is not possible to identify the signal state with 
a simple two-mode squeezed state as generated by an optical 
parametric oscillator. The remarkable feature of these states 
is that the quadrature variances are smaller than theoretically 
possible using a coherent state local oscillator of the same size 
[c.f. Eq.©], a fact only made possible by the entanglement 
between local oscillator and signal beams (for further detail 
see the Appendix.) 

Number fluctuations are an important consideration in an 
experimental setting, and will typically be at or above the 
shot-noise limit. We repeated the above simulation with an 
initial coherent state condensate with Poissonian statistics 
with a mean of 100 atoms. We do not plot the results here 
as they are essentially identical to those in Fig.|2](a-c) for an 
initial number state. Simulations of systems beginning with a 
highly-mixed initial state with number fluctuations well above 
the shot-noise limit also demonstrate entanglement. We plot 
the results for an initial thermal (chaotic) distribution with a 
mean of 100 particles in Fig. |2](d— f). We observe that in this 
case the maximal amount of violation and the range of values 
of xt where entanglement is demonstrated are both reduced 
compared to an initial Fock state, but it is present nonetheless. 
In fact, for small values of \t the violation slightly increases 
with fluctuations. This result applies to both photons and 
atoms, and may suggest that the demonstration of entangle- 
ment in a range of quantum optics experiments are relatively 
unaffected by the intensity fluctuations of the laser sources. 

In conclusion we have derived three new criteria for identi- 
fying continuous variable entanglement when local oscillators 
in arbitrary states are used in the quadrature measurements. 
We have shown that these criteria can be violated and entan- 
glement demonstrated for degenerate four-wave-mixing with 
as few as 100 particles. In this situation the direct application 
of the criteria of Duan et al. [9J], Simon C3] and Reid d 
is inappropriate, and can lead to either falsely identifying en- 
tanglement, or failing to identify entanglement when the local 
oscillator modes are sufficiently small. We have also shown 
that initial number fluctuations will not necessarily prevent the 
demonstration of entanglement in experiment. 

The authors would like to thank Margaret Reid, Ping Koy 
Lam and Claude Fabre for useful discussions, and acknowl- 
edge financial support from the Australian Research Council 
Centre of Excellence program and grant FF0458313. 



APPENDIX 

Here we expand on two issues discussed in the main text. 
Firstly we compare and contrast the new entanglement criteria 



we derive to the earlier results of Duan et al. @], Simon ifioll 
and Reid III ill when using coherent state local oscillators of 
arbitrary size, and illustrate with an example. Secondly, we 
present a more detailed analysis of degenerate four-wave mix- 
ing of a Bose-Einstein condensate as presented in the main 
text, and make a comparison for a range of particle numbers. 

The issue of measuring quadrature variances using a limited 
size but uncorrelated local oscillator was addresed by Zhang 
et al. tl5[] . In this experiment it was necessary to reduce the 
intensity of the local oscillator to be of a similar magnitude 
to the signal to avoid saturating the photodetector. As is com- 
mon in experiments in quantum optics, the local oscillator was 
separated from the beam generating the signal early in the ex- 
periment using a beam splitter. If the laser beam incident on 
the beam splitter is a coherent state (or can be written as a 
classical mixture of coherent states), then one can assume the 
signal and local oscillator beams are not correlated. We will 
examine the case where the source beam (and therefore local 
oscillator) is in a coherent state \(3). 

If the signal (local oscillator) mode has annihilation opera- 
tor a (b), the variance of the intensity difference of the homo- 
dyne measurement is 



Var[a f 6 + atf] = (tfb)Vax[a + a f ] + (a f a) 



(16) 



where we have assumed (b) = (3 is real. We can use this re- 
sult to determine the canonical quadrature from the measured 
quadrature, and therefore derive the equivalent of the Duan 
et al. inseparability criterion |§| when using small coherent 
states for the local oscillators. We find that all separable states 
obey 



Yax[X t ± X 2 ] + Var[Yi =p Y 2 ] > 



a a 



(b'A) (btb 2 ) 



(17) 



This inequality is easier to violate than the inseparability cri- 
terion we derive as Eq. (7) in the main text, but requires that 
the local oscillators be in a coherent state. 

The same corrections can be applied to derive an EPR crite- 
rion analogous to Eq. (10) in the main text, assuming coherent 
state local oscillators 



Var [A\ ± X 2 ] + Var [F x T %} > 



2 U 2/ 



(18) 

(bib,) (b\b 2 ) 

A similar analysis can be applied to the EPR criterion of 
Reid 111 ill , with final result 

'*t* v Var^,^] 2 



(ala 



2 "2/ 



AL[X 2 ]=Var[! 2 ]- 

2 

b 2 ) var[yj-fp 



(blk) Var[Xj - 
Var[F 1} y 2 



A? nf [Y 2 ] = Var[F 2 ] - , L ^ , 2j t , ; , (19) 



(btk 



(blh) 



A? nf [X 2 ] A? nf [F 2 ] >1, 
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again assuming coherent state local oscillators. 

The focus of this work is situations when it is not possible 
to assume a coherent state local oscillator of any size. The 
main text illustrated this with the example of degenerate four- 
wave mixing beginning from a number state of 100 bosons 
and using the remnants of this mode to form the local oscil- 
lators. "Perfect" local oscillators (i.e. large and coherent) for 
the quadrature measurements cannot be used to demonstrate 
entanglement in this system. This can be seen from Eq. (14) 
in the main text where the state basis clearly implies tripartite 
entanglement; using an external local oscillator is equivalent 
to tracing over mode resulting in a mixed, separable density 
matrix. Thus it is necessary to use mode as the source of our 
local oscillators in order to detect entanglement in this system, 
and the new entanglement criteria that we have derived can be 
correctly applied. 

We now illustrate the important differences between using 
perfect local oscillators and the method of quadrature mea- 
surement described in the main text with a numerical simu- 
lation. To do this we begin with a coherent state in mode 0, 
and vacuum in modes 1 and 2, represented by the state- vector 
\a, 0, 0). The state then evolves under Hamiltonian in Eq. (13) 
of the main text, and so we can write the state vector in a re- 
duced Fock basis, 



Coherent state \a = 10} 



N,m 



t (t)\N -2m,m,m). 



(20) 



The exact evolution can be computed efficiently as the Hamil- 
tonian conserves total number and so the evolution of {cN >m } 
can be calculated separately for each N. It is then a simple 
matter to extract the expectation values and variances of the 
quadratures and populations. 

To make a direct comparison with the entanglement criteria 
of of Duan et al. we rearrange our separability criterion 
(Eq. (7) in the main text) to the following form 



Var [X 1 ± X 2 ] + Var [Y\ =p Y 2 ] 



1 



1 



(6p 2 ) 



>4, 



(21) 



For a comparison with the results of Reid 111 111 , we rearrange 
the two EPR criteria (Eqs. (10) and (8) in the main text) 



V&r[X 1 ± X 2 ] + Var[Yi T Y 2 ] 



(a+a 2 ) 



> 2, 



(" 



> 1. 



(22) 



(23) 



We note that (aja^ / (6^) = {a\a 2 ) / (b 2 b 2 ) f° r our system, 
and thus that the LHS of Eqs. d2TT > and d22i i are equal. Setting 
the denominators on the LHS of Eqs. (121122123b to be one re- 



covers the inequalities in Refs. |2,LL1|]- We plot a comparison 
of the different separability and EPR inequalities in Fig. [3] 

Figure |3l a) shows the evolution of the mean population of 
the local oscillator and signal modes, which is identical for 
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FIG. 3: (Color online) (a) The population of the initial condensate 
(solid line) and the signal modes 1 and 2 (dashed line) as a function 
of time for degenerate four-wave mixing. The condensate begins in 
a coherent state with a mean of 100 particles, (b) The rescaled sepa- 
rability criterion [Eq. d21H is plotted (solid line) with a comparison 
with that of Duan et al. J9(] (dashed line), (c) The rescaled EPR crite- 
rion [Eq. d23H is plotted (solid line) with a comparison to the result 
of Reid ifliFtdashed line). 



an initial condition of either a coherent state or a Poissonian 
mixture with a mean of 100 bosons. Figures [3|b,c) show a 
comparison between using the perfect local oscillator with 
previously known entanglement criteria Jgl [111] and the cor- 
responding results derived in this work. For short times they 
are in agreement, but later the perfect local oscillator clearly 
outperforms the remnants of mode 0. However, the point of 
maximum violation is not significantly different between the 
two methods. 

We again note that for an initial Fock state, or any classi- 
cal mixture of Fock states, the perfect local oscillator will be 
unable to detect any entanglement between the signal modes. 
Also, if the state of an external local oscillator is somehow 
comprimised (e.g. by the Kerr nonlinearity, which is signifi- 
cant in BECs) then the remnants of mode may make a more 
suitable phase reference for an experiment. 

We now present a comparison of degenerate four-wave 
mixing with a BEC in an optical lattice for different initial 
particles numbers. In Fig. [4] we plot a range of results from 
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FIG. 4: (Color online) The population dynamics and entanglement 
criteria for Fock state initial conditions with 10 (solid line), 100 
(dashed line), 1000 (dotted line) and 10000 (dash-dotted line) par- 
ticles versus the scaled time N\t. (a) In each case atoms are trans- 
ferred from the condensate (blue lines) to modes 1 and 2 (red lines), 
(b) The rescaled separability criterion [Eq. d21l l. blue lines] shows 
that the amount and duration (in scaled time units) of violation in- 
creases with particle number. Furthermore, the violation reoccurs af- 
ter the local oscillator becomes smaller than modes 1 and 2. The un- 
depleted pump solution is plotted for comparison (red line), (c) The 
rescaled EPR criterion [Eq. {23}, blue lines] shows similar behavior. 
The undepleted pump solution is also plotted (red line). 

N = 10 to TV = 10000 particles. For a clear comparison we 
have plotted the results on a time scale N\t for which the ab- 
solute populations of modes 1 and 2 are similar for all intial 
atom numbers at small times. In this region the population of 
mode is approximately constant with a 2 , as aj 2 w N, and so 
the Hamiltonian is 

H w i\N (a^ - a\aU . (24) 

This is the undepleted pump approximation with the exact so- 
lution in the Heisenberg picture 

a^t) = coshiNxtfa^+sinhiNxtfaliO), (25) 
a 2 {t) = cosh(N X t)a 2 (0) +smh(N X t)a\(0). (26) 

This allows us to compare with the limit TV — > oo; we can also 
assume that mode will act as a perfect local oscillator in this 



limit and apply the entanglement criteria of Duan et al. |@1 and 
Reid fllll . 

Figure Ufa) shows that the relative number of atoms trans- 
ferred from mode to modes 1 and 2 remains insignificant 
for a longer timescale N\t for larger initial populations. In 
Fig. 2|b) and (c) we plot the rescaled entanglement criteria 
[Eqs. ( 1211221231 )1, and observe that the entanglement is de- 
tected for a larger range of N\t for larger values of N. For 
comparison we have plotted the original entanglement crite- 
ria of Duan et al. @| and Reid fllll for the undepleted pump 
solution in red. The finite particle number solutions begin to 
deviate from this curve around the same time that the popula- 
tion of mode decreases by a few percent. 

In Fig. @Jb) the separability criteria is violated in a region 
where modes 1 and 2 are larger than the local oscillator mode 
for N > 100 particles. This violation increases with in- 
creased particle number, but we do not demonstrate the EPR 
paradox even with N = 10 4 particles (corresponding to our 
computational limit). The state of the system at this time is 
complex and dissimilar to the two-mode squeezed state pre- 
dicted by the undepleted pump approximation. It is interest- 
ing to note that at these times the system exhibits quadrature 
variances that are much smaller than is theoretically possible 
using a coherent local oscillator of the same size. We can see 
from Eq. ( [Tol l that using a coherent state local oscillator sets a 
lower bound to the measured quadrature variance. This is also 
true for the sum or difference of quadratures 

Var[X 1± ! 2 ]>M^ + MR (27) 
J " (bib,) (b\b 2 ) 

The variances of the quadrature sums or differences predicted 
in our simulations are as small as half of this limit. The entan- 
glement between the signal and local oscillator modes allows 
a reduced level of uncorrelated fluctuations. 

Some time after the recurrence period the entanglement in- 
equalities are again no longer violated once the population 
transfer from mode to modes 1 and 2 has begun to reverse. 
At this point the quantum state of the system has evolved quite 
significantly, and we observe interference fringes in the num- 
ber probability distribution for each of the modes. These con- 
tinue to become more complicated, and entanglement will not 
be observed at any later time. 
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